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ABSTRACT 


A new method is developed by which certain functions of the noncommuting operators 
A and B can be expressed in a systematic and efficient way in terms of successively 
higher commutators of the operators. The procedure is designed primarily for poly- 
nomials, exponential functions, and analytic functions of the operator (A + B). The 
method is related to the ordered-operator calculus of Feynman. Certain commutation 
properties of the operators A and B are represented by polynomials in A and B called 
iterated commutators. The functions are converted into equivalent forms involving the 
iterated commutators. All iterated commutators of degree (j + 1) in A and B are com- 
bined in a suitable sum to form a new operator G^. By working only with the Gq^ and 
their commutators, the calculations are simplified and the results are given in an ex- 
tremely compact notation. 


ITERATED COMMUTATORS AND FUNCTIONS OF OPERATORS 
by Howard C. Volkin 
Lewis Research Center 

SUMMARY 

Certain commutation properties of two operators A and B are represented by poly- 
nomials in A and B called iterated commutators. A method is developed for system- 
atically expressing functions of A and B in terms of the iterated commutators. The 
method is related to the ordered-operator calculus of Feynman. Compact expressions 
result because all iterated commutators of the same degree always appear in the same 
linear combination. The procedure is designed primarily for polynomials, exponential 
functions, and analytic functions of the operator (A + B). Applications to the important 
case of exponential functions are given. 


INTRODUCTION 

Functions of noncommuting operators have been studied extensively over the years, 
usually within the context of problems arising in mathematics and quantum physics. A 
survey of results, along with numerous references, has been given recently by Wilcox 
(ref. 1). One approach in working with such a function is to express it as an expansion in 
successively higher commutators of the operators involved. The purpose of this report 
is to present a new method by which explicit forms of the desired type can be developed in 
a systematic and efficient way for certain functions of interest. 

This method is related to the ordered-operator calculus of Feynman (ref. 2). In 
Feyman’s operator calculus, an ordering parameter indicates the sequence in which non- 
commuting quantities appearing an any product are to be taken. An operator A, which 

does not otherwise depend on the ordering parameter s, may be written A , with the 

s 

ordering index attached as a subscript. Then it is understood that the form of A does 

s 

not depend on s, but its order of operation does and is defined by the convention: The 
operator with smallest index operates first, and so forth. Thus, A B. equals BA when 
s’ > s, equals AB when s > s', and is undefined when s = s'. In the more general 
case of an operator that is a function of the time, the convention is that time is the 


ordering parameter. We can use the notation A(s) to include either case. 

With the ordering parameter convention, any functional of the operator functions 
A(s), B(s), . . . has a unique interpretation, because the order of operations in any situ- 
ation is automatically specified by the parameter values. In manipulating such a func- 
tional, the algebra and calculus of ordinary numbers may be applied, as though A(s), 
B(s), . . . were ordinary commuting functions. "Disentangling" is the name given by 
Feynman to the process of finding an operator expression for the functional in the usual 
positional notation, where the order of operation is represented by the position in which 
the operators are written. Simple examples of disentangling are the following: 


exp(Aj + Bq) = exp(A) exp(B) (1) 



An example of a full functional is the well-known Dyson formula for the Neumann- 
Liouville expansion of the unitary time development transformation. The transformation 
is an exponential functional of a time -dependent Hamiltonian operator function. The 
ordering is usually exhibited explicitly by the chronological operator, which may be 
omitted when the Feynman notation is adopted. 

The Feynman formalism permits all the results of ordinary analysis to be used in 
handling operator functionals. However, no general procedure for disentangling func- 
tionals is available. Moreover, there does not appear to be any direct way of exploiting 
any simple commutation properties that the operators may happen to possess. The 
method presented here deals with this second aspect. Certain commutation properties of 
two operators A and B are represented by polynomials in A and B called iterated 
commutators. The procedure given converts polynomials, exponential functions, and 
analytic functions of (A + B) into equivalent forms in terms of the iterated commutators. 
All iterated commutators of degree (j + 1) in A and B are combined in a suitable sum 
to form a new operator G^. By working only with the and their commutators, the 

calculations are simplified and the results can be given in an extremely compact notation. 

In the next section, the iterated commutators are defined and their relevant prop- 
erties are given. In the section EXPANSION IN ITERATED COMMUTATORS, relations 
are developed that convert a suitable function of two operators into the desired form. 
Applications to the important case of exponential operators are given in the section 
EXPONENTIAL FUNCTIONS. As the derivations required are generally quite straight- 
forward, complete details are given only where deemed necessary and a suitable method 
is indicated in all other cases. 
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SYMBOLS 


A 

A s 

A(s) 

B 

C 

C (i) 

E (i) 

F(x) 
G/ 


'(j) 


H 


H (i) 

H„ 


p 0) 

R 


operator 

operator with ordering index s 

operator which may depend on ordering parameter s 
operator 

operator in eq. (29) 
operator in Zassenhaus formula 
operator in expansion of C 
analytic function of x 

linear combination of iterated commutators 
operator in BCH formula 
operator in expansion of H 
A + B 

F[MA + B)] -1 F[X(A 1 + B q )] 
operator in expansion of P 
F[A(A + B)] - F^Aj +■ B q )] 


s ordering parameter 

{(A) n ,(B)} [A, . . . [A,[A,B]]. . . ], n commutations 

{(A),(B) n } [. . . [[A,B],B]. . . ,B], n commutations 

{ (A) m , (B) n } L { ( { (A) m , (B) } ), (B) n “ 1 } 

{ (A) m , (B) n } R {(A)" 1 " 1 , ({ (A), (B) n } ) } 


Superscript: 

t 


Hermitian adjoint 


ITERATED COMMUTATORS 

The simple iterated commutator is 

{ (A) n , (B) } s [a, ... [A, [A,B]]. . .] 


(3) 
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which defines a linear operation on B obtained by n successive commutations with the 
operator A. The form of equation (3) is a polynomial whose terms each have n factors 
of A and one of B. Similarly, 

{ (A), (B) n } s[. . .[[A,B], B]. . B] (4) 

The case n = 0 is specified by 

{ (A)°, (B) } = B, {(A), (B)°}=A 

It is seen that 

{ (B) n , (A) } = (- l) n { (A), (B) n } (5) 

Clearly, the iterated commutators (eqs. (3) and (4)) arise when two operators, one 
of which is raised to a power, are transposed. The following four theorems are estab- 
lished by induction: 


{ (A) n , (B) } = 2 (-l) r (J)A n - r BA r 
r=0 

[A n , b] = £(J) « A ) r , <B)}A n - r 
r=l 

[A,B”]= X;(J)B n - r {(A),(B) r ) 
r=l 

n 

{ (A) n , (BC ) } = ^ (r)t ( A ) n_r » ( (A) r , (C) } 

r=0 


( 6 ) 

(7a) 

(7b) 

(8) 


A familiar example in which these iterated commutators occur is the Lie expansion, 
which is the power series 

exp (A. A) B exp(-XA) = (X r / r!) { (A) r , (B) } 

r=0 
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I 


There is no unique way to prescribe a commutator {(A) m , (B) n } that reduces to 
the forms in equations (3) and (4) for the special cases n = 1 and m = 1, respectively. 
Various generalizations suggest themselves. The extended definition of iterated commu- 
ator is chosen to be 


{ (A) m , (B) n } L = {( { (A) m , (B) } ), (B f- 1 } (9) 

hi the "L-type” of extended commutator (eq. (9)), the simple commutator (eq. (3)) with 
the m-fold iteration of the operator A (the operator appearing in the left-hand position) 
is formed first, followed by the remaining (n-1) iterations of the other operator, B. An 
R-type commutator, for which the prescription is reversed, is easily seen from equa- 
tion (5) to satisfy the relation 

{ (A) m , (B) n } R = (-l) m+n-1 { (B) n , (A) m } L (10) 


The iterated commutators (eqs. (9) and (10)) consist of a sum of products each having m 
factors of A and n factors of B; they describe the commutation of a power of one 
operator with a power of another operator. From equations (6), (7), and (8) follow 
directly the reduction formulas 


{ (A) m , (B n ) } 
{ (A n ), (B) m } 


n m 

£ £ (;)(T:i 1 )f<A) m - t ,(B''- r )K(A)',(B) r ) R 

r=l t=l 

n m 

E E (?)( r ?:i 1 )HA) r ,(B) t } L {(A n - r uB) m - t } 

r=l t=l 


( 11 ) 

( 12 ) 


Equations (7a) and (12), for example, convert any polynomial in A and B into a sum 
of terms, each of which is a product of iterated commutators and an ordered product 
A m B n . 

If A is the differentiation operator with respect to a parameter on which the oper- 
ator B depends, then equation (3) becomes the n th derivative of B, that is, 


d n B 

dt n 



(13) 
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and equation (9) describes the commutation properties of the derivatives with the powers of B, 



(14) 


Employing in turn the identities [A, B]^ = -[A^, B^] and equation (10) gives the following 
for the Hermitian adjoint: 


{ ( A) m , (B ) n } j= (-lr 4 *- 1 { (Atr, (B t) n } L 

{(A) m ,(B) n } L t={(Bt) n ,(At) m } R 


(15) 

(16) 


EXPANSION IN ITERATED COMMUTATORS 

First some useful disentangling expressions involving the quantities (A^ + Bq) 11 will 
be developed. These "ordered powers" of (A + B) are simply 

(Aj ♦ B / = £ (J) A n-r B r (17) 

r=0 


that is, they are obtained as ordinary binomial expansions in which A^ and Bq are 
treated as though they commuted. When each term is written in the usual positional con- 
vention, powers of A appear always to the left of powers of B. Multiplying both sides 
of equation (17) from the left by (A + B) gives 


(A 2 + B 2 )(A 1 + B 0 ) n = (Aj_ + B 0 ) n+1 


n 


-E (,? 1 )[*" 1 ' r . B > r - 1 


(18) 


where, for convenience, both notations have been used on the right. After using, in suc- 
cession, equations (7b), (6), (7a), and (9), the second term on the right-hand side of 
equation (18) becomes 

n n+l-r r-1 

£ E E <- 1 > t (z--i)( n+ t' I ')( r p 1 ) An+ 1 - r - t B r - 1 -e t <A)‘,( B r 1 } L 

r=l t=l p=0 
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The substitutions j = p + t, k = t - 1 and l =n+l-r-t separate the summation into 
independent sums over i and k. The result allows equation (18) to be written in the form 

n 

<A 2 + B 2 )(A l + B 0 ) n = (Aj + B 0 ) n+1 - £ (") (Aj + B/-i G fl) (19) 

j=l 


where 


j-1 




k=0 


( 20 ) 


The operator is a homogeneous polynomial of degree (j + 1) in the operators A 

and B. The subscript j is bracketed to signify that it has no ordering significance. It 
is found, for example, that 


G(1) = [A, B ] 

G (2) = 2 { (A), (B) 2 } - { (A) 2 , (B) } 

G (3) = 3 { (A), (B) 3 } - 3 { (A) 2 , (B) 2 } L + { (A) 3 , (B) } 

G (4) = 4{ (A), (B) 4 } - 6 { (A) 2 , (B ) 3 } L + 4 { (A) 3 , (B) 2 } L - {(A) 4 ,(B)} 

On the right-hand side of equation (19), all the iterated commutators (eq. (9)) whose 
indices m and n add up to the same value appear in the single operator (eq. (20)) for 
which j = m + n - 1 and in none of the others. If only a finite number of iterated com- 
mutators are nonvanishing, then all G^ with values of j larger than a certain value 
will vanish. As an illustration, consider the case where A is proportional to p 2 , B is 
a polynomial in x of degree l , and [p,x] = c, an ordinary number. With the help of the 
identity 


m 

[VV] = £ (“) i a - D- • • u - r ♦ i)cV - r P m - r 

r=l 

it is found that 
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o rn 7 m 

{ (p 2 ) ,(**)} = Z q r ( x )p 

r=0 


where q r (x) is a polynomial in x of degree l - 2m + r. The iterated ^commutator 
vanishes if m>Z. On the other hand, when m { (p^) m , (x^)"* m+ }l vanishes 
if j — 2m + 1. Hence, when j & 21 + 1, each left commutator (eq. (9)) appearing in 
vanishes and = 0. 

With the help of the recursion relation (eq. (19)), the following disentangling formula 
can be verified directly by induction: 

(A x + B 0 ) n - (A + B) n 


n-1 

-E 

r=l 


+i)(A + B) n-1-r G 


n-2 


(r+1 


(r) 


+ Z E 

r=2 i+j=r 




(i) G (i) 


+ 


11-0 

E E 


r=3 h+i+j=r 


(r: 3 )( r J +2 )( h T 1 ) <A + B,n ' 3 ' r ' 


(h) G (i) G (j) 


n-4 




r=4 g+h+i+j=r 


+ S 


N 


( 21 ) 


The number N of summation terms on the right equals the largest integer that is less 
than or equal to (l/2)n. The final term Sjg- is (n - l)(n - 3) . . . (5)(3) [ G (1 ,] <1/2)n > 
when n is even, and is the sum term for r equal to (1/2) (n ± 1), when n is odd. A 
dual form of equation (21) holds, with ordered powers of (A + B) appearing on the right 
in place of ordinary powers: 
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(Aj + B 0 ) n - (A + B) n 


S (r+l)( A l + B (P G (r) XI XI (jr+zX**?" 1 ) ^ A 1 + B 0^ 2 ^(j ) G (i) 

r=l r=2 i+j=r 


n-3 


E E ( r ;3X r rTi +1 )( A i +B o)“- 3 - r G 0 )Oftfo.) + --- + <- i ) N - ls N w 

r=3 h+i+j=r 


The form of the dual expression (eq. (22)) is similar to equation (21) except that (1) the 
order of the subscripts on the factors is reversed and (2) the sign on successive 

sum terms alternates, beginning with a plus sign on the first term. 

The order of the operator factors in each term of equation (21) can be reversed by 
taking the Hermitian adjoint of each side. The Hermitian adjoint of a functional that is 
written in the Feynman notation is gotten by replacing every operator A(s) by A^(s) and 
exact ly reversing the ordering convention on the parameters of the adjoint operators. 

For a simple function of operators, the result is clearly F(A^, Bq)^ s= F(A^, b|). Equa- 
tion (15) shows that is obtained from the form (20) by substituting At,Bt for A, B 

respectively, and inserting an overall factor of (-1)- 1 . If A and B are Hermitian, 
then = (-l^G^. 

The application of equations (21) and (22) to analytic functions of (A + B) is straight- 
forward. Let the analytic function be 


F(x)= X]4 f n Xn (23) 

n=0 n! 

In the formula 

f[\(A + B)] = f[\(A x + B 0 )] + R (24) 


the term R has the series expansion 

R = (sr)^ 1 * + (^) r < 2 > * (jt) r < 3 > + " • 


(25) 


9 



Equation (22) gives the operator directly in terms of the constant f n , ordered 

powers of (A + B), and the G^. An alternative formulation is the disentangling formula 

f[x(A x + B q )] = f[x(A + B)]P (26) 


where 


P = 1 + 


(§ 




(i) 


( 3 ) 


(27) 


Equating the coefficients of like powers of X in equation (26), with the help of equa- 
tion (21), yields the operators P^. The first few terms are 

P (2) = (^) G (2) + - f l f 2>< A + B >°(1) 

P (s) ■= ^G (s) + 4 - fjfjX A + BX5 (2) + M5 S [«*f 4 - ipf 2 - 2f ^3 - fjf^fA + B)^ 


A dual form of equation (26) can also be developed, namely, 

f[x(A + B)] = F[x(A 1 + B 0 )]P _1 (28) 

where P” 1 is expressed in terms of ordered powers of (A + B). Other expansion 
methods for analytic functions of (A + B) and references to earlier work in this area 
are given by Kumar (ref. 3). 


EXPONENTIAL FUNCTIONS 

The operators P^ in equation (27) take a particularly simple form when F is the 
exponential function, in which case f n = 1 for all values of n. This leads to a compact 
statement of the theorem 


exp(XA)exp(XB) = expX(A + B)exp(C) 


(29) 
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where 


K*) e < 1 > + (*) e < 2 > + (*) 


E (3) + • • • 


(30) 


If equation (29) is written in the form 

exp [MAj + B 0 )] - exp [x(A + B)] = exp [x(A + B)J[exp(C) - l] 


(31) 


the coefficient of (X n /n! ) in the power series expansion of the left-hand side is given by 
equation (21). The corresponding coefficient on the right is 


N n-m 

EE E 

m=l r=m a+b+. . . +i+j=r 


1 ( n \/r+m\ /h+i+j+3Vi+j+2\ 

~ rVr+mAa+1/- ' A h+1 Ai+1 ) 


x (A + B) n - 1 --”E (a)E(b) . . . E(h) E (i)E(i) 


(32) 


E/an - + [p(i)>^(2)] 


In the m product term of the previous sum (eq. (32)), there appear m factors of 
and m binomial coefficients. The are found by equating coefficients of the same 

power of X on each side of equation (31). The results through ninth order in X are 

E (i) = G (i)’ i = 1 ’ 2 ’ 3 

: (4) = °(4) 

E (5) = G (5) + 0[ G (1)’ G (3)] 

E (6) = G (6) + (f)[ G (D’ G (4)] + (f)[ G (2)> G <3)] - (£) {( G d)) 2 ’ G (2)} 

E (7) = G (7) + 7 [ G (1)’ G (5)] + 7 [ G (2)’ G (4)] + (y){( G (1))’ ( G (2)) 2 } 

E (8) = G (8) + 10 [ G (l)’ G (6ll + 14 [?(2)’ G (5)] + 7 [ G (3)’ G (4)] “ 12 {( G (1)) 3 ’ (°(2))} 

+ 3 {( G (1))M G (4))} + 18 [[ G (1)’ G (2)]’ G (3)] + 10 [[ G (1)’ G (3]J G (2)] 
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Note that each of the terms occurring in is essentially predetermined to 

within a numerical coefficient. Consider any single term in the sum of terms that 
makes up the complete expression for a given E^. If the subscript on every 
appearing in the term has its numerical value increased by unity, then the numbers 
so obtained add up to (j + 1), the degree of the iterated commutators in y 
The formula of Zassenhaus (ref. 4) 


exp|x(A + B)j = exp (a A)exp (xB )exp ( ©Jexpf — 


2 ! 


3! 


( 33 ) 


corresponds to the complete factorization of exp(-C) in the inverse form of equation (29). 
With the help of equation (22), the operators are obtained by equating coefficients 

of like powers of x. Again, the results out to the ninth order in X are 


C (D = - g (d> 


i = 1,2,3 
C (4) = - G (4) - 6 [ G (1)’ G (2)] 


C (5) “ -G (5) ' 10 |p(l)" G (3)] 


C (6) = - g ( 6> - 15 [ G (1)- G (4l] - 20 [ G {2)- G (3>] - 45 (( G l) 2 > < G 2’} 

C (7) = - G (7) - 2l [ G (l)> G (5)] - 35 [ G (2V G (4)] - 210 {( G (1)> ( G (2)) 2 } - 105 {( G <1>) 2 ' ( G (3))} 
C (8) = - G (8) - 28 [ G (1)> G (6)] - 66 [ G (2)> G (5)] - 70 [ G (3)- G (4)] ' 210 {( G (1)) 2 > ( G (4))} 

-420{( G (l)) 3 , ( g ( 2))} + 420 [|p(l)’ G (2)]> °(3)] + 560 [l G (l)’ G (3)]' G (2)] 

Because of the Jacobi and other identities, the expressions given for the operators 
E^) and are not altogether unique. However, they display concisely the iterated- 

commutator nature of the terms in the expansions of equations (30) and (33). No attempt 
has been made to obtain general expressions. Nevertheless, for certain types of terms, 
the general form of the numerical coefficients can be inferred from the specific values 
already given. For example, contains the term (-l/2)j(j - 1 )Jg^^,G^_2^ . 

This procedure may also be applied to the Baker-Campbell-Hausdorff (BCH) formula 


exp(xA)exp(XB) = exp(H), 


(34) 
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where 



In the case of the BCH formula, the method offers no great advantage over previous 
methods (refs. 4 and 5) that have been developed specifically for this problem. However, 
the can be computed within the framework of the procedure given herein without 

introducing additional elements (such as operator derivatives or auxiliary operator func- 
tions) that the more specialized methods employ. Equating the coefficient of (A n /nl) in 
the expansions of fexp(AH) - exp(AH Q )] to the expression given by equation (21) yields 

H (i)= G d) 

H (2) = G (2) + (|)[ H o’ G (1)] 

H(3) = G (3) - 2 [ H o> G ( 2 {] - {( H o) 2 ,(G(1))} 

H <4> * G (4) ^ (|)[ H o. G (3)] - [ G (X)> G (2)] ^ (|){(H 0 ) 2 , ( g (2) )} * (f){(H 0 ), (G (1) ) 2 } 

H (5) - G (5) * 3[H 0 ,G (4) ] ♦ (!)[g (1) ,G (3 ^] ♦ (|){(H 0 ) 2 ,(G (3) )) + ^,0^, G (2) ] 

4 2[[H„. G ( 2 )], G ( 1)] 4 {(H 0 ) 2 ,(G (1) ) 2 } L - (i){(H 0 )M G (X))} 

When expressed in terms of the G^, the contain considerably more terms 

than either the or the The additional terms are the ones containing H Q . 

Moreover, the relations^ 

G (2) = {< A >. < B > 2 } 4 [o (1) , H 0 ] 

and 

■^No such simple relations exist for higher Gq), however. 
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0(3) = {(A), <B) 3 } + [g ( 2) , H d ] 


can be used to reduce the expressions for the H^. For example, the following can be 
written: 

H(2) = (|){(A). (B) 2 } - (|)(( A ) 2 . <B)J 

“(3) = {< a > 2 > < b > 2 }l 

There exist any number of intermediate formulas between the BCH form (eq. (34)) 
and equation (29), all of which can be written as 

exp(xA)exp(xB) = exp XH q +^ H ( 1) + • • • + ^ H (n-1) ex P n) + ’ ' * (36) 

with equation (29) by n = 1 and equation (33) by the limit n — The operators on the 
right-hand side of equation (35) can be evaluated by the procedure given herein. In 
particular, the operators in the first exponential on the right are those of the BCH 
formula, as follows directly from the fact that the right-hand-side expansions of equa- 
tions (34) and (36) are the same up to the n^ power of X. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, August 21, 1968, 

129-02-07-07-22. 
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Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and Notes, 
and Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION DIVISION 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 



